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Abstract 

The Hartree-Fock-RPA approach is applied to the ID anti-ferromagnetic Heisen- 
berg model in the Jordan- Wigner representation. Somewhat contrary to expecta- 
tion, this leads to reasonable results for spectral functions and sum rules in the 
symmetry unbroken phase. In a preliminary application of Self-Consistent RPA to 
finite size chains strongly improved results are obtained. 
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The ID Anti-Ferromagnetic Heisenberg Model (AFHM) [T] belongs to one of the most 
classic many body research fields. It was the first many-body model to be solved exactly 
by Bethe with his famous Bethe ansatz solution [2] and it has been a playground for 
testing and developing many body approaches ever since. In spite of tremendous progress 
in the understanding of many facets of the model, it continues to be a very active field of 
interest and research [3J H] . 

In an attempt to apply extended RPA-like theories, as e.g. the recently developed 
Self-Consistent RPA (SCRPA) approach j5j[6], we found out, somewhat to our surprise, 
that even the standard RPA theory has so far not fully been developed in its application 
to the 1D-AFHM. This probably stems from the fact that the HF-RPA scheme is usually 
considered as unreliable in low dimensions [7] . However, it will be found that RPA shows 
interesting features even in this ID-version of the model. On the other hand SCRPA 
shows very promising results in a number of cases [5] and in a first application for a small 
numbers of sites we here get very close agreement with results from exact diagonalization. 

We mostly consider the isotropic Heisenberg model with anti-ferromagnetic coupling, 
however some consideration will also be given to the anisotropic case, for instance to the 
XXZ-anisotropic Heisenberg Hamiltonian 

where S l n is the spin operator of the n-th site, g is the anisotropy parameter and N is the 
number of sites. Periodic boundary conditions (S l N+1 = S[) are applied. We will work 
with the AFH model in the Jordan- Wigner (JW) representation [8]. The Hamiltonian in 
momentum representation is given by [9] 

H = ^ + J2 £ °k^ k+ \ Yl ^^ksk^l^i^tpks (2) 

k k\k2k^ki 
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with 

e° k = cos(k) - g (3) 

and 

v kl k 2 k s k 4 = j-f (cos(&i - k 3 ) - cos(k x - fc 4 )) + k 2 -k 3 - k A ). (4) 
The isotropic point is given by g = 1 whereas for g ^ 1 we have the XXZ model. In order 

7T 

to take care of the boundary conditions the momentum sums run over the values k = j— 
with : 

j = ±1, ±3, ±5, . . . , ±(N — 1) for n even 

j = 0,±2,±4,...,±(JV-2),JV for n odd (5) 

where n is the fermion number (in this note we only consider the sector n = N/2) and 
where we have assumed, without loss of generality, that N is even. The Kronecker symbol 
is defined by : 

5(p-p) = 1 if P = P '±27tt; t = 0,1,2,- ■■ 

= otherwise. (6) 

We first want to study the Hartree-Fock (HF) approximation, keeping translational 
invariance, i.e. momentum as a good quantum number. The JW-HF single particle 
energies are given by e k = [H, ^/] j). Here [.., ..] stands for commutator, {.., ..} for 

anticommutator, and (...) denotes the expectation value in the HF state. One obtains 

ek = £ * + % (1 - cos(A; - />)) . (7) 

h 

where the sum goes over the occupied states "h" only. The HF groundstate energy is then 
given by 



< F = ^+>:i4+^>:(i-cos(/,-/,'))i. (8) 

It is interesting to note that with respect to the Neel state the HF energy of the JW 
transformed isotropic AF Hamiltonian has actually a lower minimum. In table Q] we give 
the HF energy for the Neel state and the JW-HF energies as a function of the number of 
sites iV for the isotropic situation with g = 1. 



h \ h> J 



N 


4 


6 


8 


12 


Neel energy 


-1.0 


-1.5 


-2.0 


-3.0 


JW-HF energy 


-1.9142 


-2.6667 


-3.4667 


-5.1077 


Deformed JW-HF energy 


-1.9142 


-2.6667 


-3.4856 


-5.1927 



Table 1: Mean- field energies with various approximations for some low number of site 
cases (see text for details, e.g. 'deformed JW-HF energy' stands for JW-HF energy in the 
symmetry broken phase). 
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So the JW-HF theory takes into account more of the interaction energy than the Neel 
state in real space. It seems, however, not easy to analyze the physical content in real 
space of the JW-HF groundstate. 

We now go one step further and investigate the Random Phase Approximation to 
describe excitations on top of the before considered JW-HF groundstate. For this we 
elaborate the RPA equations via the equation of motion (EOM) method [TU1 [TT]. We 
introduce an RPA particle (p)-hole(h) excitation operator 

ph 

with \v) = Ql\0) and \u), |0) the excited state and groundstate respectively. With the 
vacuum condition Q u \0) = 0, one arrives as usual [12] to the RPA eigenvalue equation 

(10) 

with 

Aphjh' = ( [ii^P > \ H > ^t' ] ) 

= (tp — £h)8 P p'5hh' + v p h>,hp' (11a) 
B ph ,p>h> = -(bP{4>p,l H >4 , h>i J p'}}) = Vpp'M' 

We get the standard RPA from Eqs. fllOp in evaluating expectation values with the before 
elaborated JW-HF groundstate leading to the expressions given in Eq. liTTl) . We show 
in Fig. CD for N = 180 and g = 1 the RPA eigenvalue spectrum as a function of the 
momentum transfer |g|. We clearly see that there is a lower branch slightly detached from 
the more or less dense remainder i.e. the p-h continuum. Actually iV = 180 is already 
quite close to the thermodynamic limit which we also calculated using the Green function 
method following closely the work of Todani and Kawasaki [T3j . Here we only show the 
result in Fig. [2j We again see that a lower energy state (lowest broken line) becomes 
detached from the continuum, the latter being embraced by the two upper broken lines. 
The continuum was also studied by Todani and Kawasaki [13] but for some reason they 
did not discuss the low-lying discrete state which is an important and prominent feature. 
In the thermodynamic limit the dispersion equation of the low-lying branch can be found 
from (see Ref. [13] Eq.(34)) 

l-f(q,u)ft(G°(q,u)) = (12) 
where G°(q,u) is the renormalized particle-hole propagator 
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G\q, w) = n — n 1^ I n — * . I » ( 13 ) 



p7TSin(|)2x/r^2 \ y/T^X^ — 



sin 



and 

f(q,u) = 2cos(g) - 2px 2 (14) 
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Figure 1: RPA Excitation spectrum for iV = 180 and g — 1. 

with x = — and p = 1 H — . From Eq. (Tl2l) also the lower u>l and upper uju 

2psm(|) 7r 

boundaries of the continuum can be obtained 

2 

uj l = (1 + -) sin(g) (15a) 

7T 

^ = 2(l + -)sin(f). (156) 

In the exact two-spinon case there is no discrete state and only a continuum limited by the 

upper (uju = 7rsin(|)j and lower (ul = — sin(g)j full lines of Fig. [2] exists [TU 0]. The 

exact two-spinon dynamic structure factor diverges at the lower limit [TU H] . Nevertheless, 
globally, the RPA spectrum for the infinite chain resembles the exact solution, even though 
locally, especially at the lower edge, things go wrong even qualitatively. For instance, 
the fact that the excitation spectrum exhibits a low lying discrete state detached from 
the continuum is of course an artefact of the RPA approach. This fact has also very 
briefly been mentioned in Ref. !T5] but no details are given there. We want to point out, 
however, and this is the main message of the present paper, that the RPA spectrum is 
already qualitatively quite similar to the exact one as seen in Fig. [3] where we show the 

7T 37T 

structure functions for q = — and q = — corresponding to RPA and exact two-spinon 

calculation [1]. The detachment of the discrete state is in fact only very mild and we 
therefore conclude that RPA may be a good starting point for more elaborate theories, 
even for situations where, a priori, it should not work well as for instance in ID situation, 




Figure 2: Excitation spectrum for the infinite chain. The lowest broken line corresponds 
to the discrete RPA branch. The two upper broken lines limit the RPA p-h continuum. 
The two full lines limit the continuum in the exact two-spinon case [I] where there is no 
discrete state. 



as considered here. 

Let us also mention that in the attractive (ferromagnetic) case quite naturally a dis- 
crete state, detached from the continuum, exists, even in the exact solution [T6l [15] . 

To understand the reason why there exists in RPA a discrete collective state detached 
from the continuum it is instructive to solve Eq. ffT2"j) for the dispersion relations cu(q) 
in the case of a finite number of sites. In the discrete case f(q,uj) actually has a more 
complicated structure and it only reduces to Eq. ( fT^l) in the thermodynamic limit. In 
order not to overload the figure we choose a somewhat lower value for N than in Fig. [TJ 
We see in Fig. H]the graphical solution of Eq. ( jT2l) for iV = 80. The structure is very similar 
to other schematic RPA models (see e.g. |llj) with the only difference that the smooth 
effective interaction is momentum and energy dependent what stems from the fact that 
here we have a 3-rank separable interaction instead of the usual one rank separability. The 
collective state gets detached from the more or less unshifted ph-roots precisely because 
it is situated at the edge of the continuum. One also can make a guess what may happen 
in the exact case. Due to screening the corresponding effective interaction f(q,u) will 
bend much more strongly at the lower edge of the ph-continuum in such a way that the 
collective state sticks closely to the last lower vertical asymptote. 

The similarity of the RPA structure function with the exact one can also be inferred 
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Figure 3: Exact two-spinon and RPA dynamic spin structure factor, in the thermodynamic 

limit, for q = — and q — — . 
H 2 H 4 



from studying the corresponding energy weighted sum rules. Let us define 

1 f 00 IE 

— uS zz {q,uj)duj = — — ^(l-cos(g)) (16a) 
1 f°° c 

7T uSi 2 J(q,u)du} = — k (1 -cos(cj)) (166) 
2tt Jo 71 
1 Z" 00 1 

— / uS^ A (q : u)dw = -(1 - 008(g)) (16c) 

27T Jo 7T 

with C = 0.72221..., k = 0.9163... and = — iV(log2 — |) being the exact ground-state 

(2) 

energy. S zz (q, u>) is the exact dynamic spin structure factor, S zz (q, u) is the exact 2-spinon 
part of the exact dynamic spin structure factor [1], and S RPA (q,u) is the corresponding 

7T 

structure function in the present RPA approach. For q = — , in RPA, the sum rule 

gives 0.318309 and the contribution to the sum rule from the discrete state is 0.226843 

3vr 

(71%) and from the continuum it is 0.091467 (29%). For q = — , in RPA, the sum rule 

gives 0.543389 and the contribution to the sum rule from the discrete state is 0.244919 
(45%) and from the continuum it is 0.298470 (55%). From the above we see that the RPA 
energy weighted sum rule Eq. f ll6cl) within ~ 50% of the exact two-spinon result Eq. ( \16ch 
and within ~ 8% of the exact result Eq. (116fej) . This latter result again underlines the 
semi-quantitative correct behavior of the RPA. 
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Figure 4: Graphical solution of the dispersion relation for N = 80 and q = — . 
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On Figs. [Hand [2] we also see that RPA actually becomes unstable for q > q c ~ 0.917T 
because the discrete branch touches zero at q ~ q c . An RPA eigenvalue coming down to 
zero inevitably means that the corresponding HF stability matrix [111 H2] 

SM=(i B A ) (17) 

also has a zero eigenvalue at the critical value q = q c and negative eigenvalues for q > q c , 
clearly indicating that the translationally invariant JW-HF solution is unstable. 

The reason why the RPA performs relatively well is probably precisely because we 
forced the system to remain in the symmetry unbroken, i.e. 'spherical', phase as it is 
well known that in ID no spontaneously symmetry broken phase exists. None the less, 
we were curious to see what HF and HF-RPA gives in the symmetry broken phase. Not 
unexpectedly we will find that deformed HF still lowers the energy (see table [T]) but HF- 
RPA leads to an unphysical excitation spectrum. Follow some details of our procedure. 

We seek a new HF basis with a stable minimum in introducing new quasiparticle 
operators which are a superposition of particle (p) and hole (h) operators in the old basis, 
that is [I7j 

Oih = UftVft - Vfc^jUr ( 18 ) 

with u 2 h + vf t = 1. Our ansatz Eq. flTSl) is motivated by the fact that the mode which 
becomes unstable occurs at q = it and therefore, suggests the order parameter {ijr h+ .^h) 7^ 
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0, involving breaking of translational invariance. Following the standard procedure for 
obtaining the amplitudes Uh, Vh in Eq. (JT%|) we express the HF-energy in the new deformed 
state 



E 



HF 




gN 
4 



+ E [ 6 h + -j^ V h,h',h,h' ) 
h \ h' ) 



Vh+Tr,h',h+TT,h' — Vh,h',h,h> 



h h 
^ ' - 2 I 2 



+ - (Vh+TT,h'+TT,h+Tr,h'+Tr + Vh,h',h,h' ~ h+TT,h' ,/l+TT ,h' ) v h'j v h 
+ '^^(Vh+w,h'+n,h,h> + Vh+TT,h>,h,h'+n) UhVhUh'Vh' (19) 



hh' 

and minimize with respect to Uh, Vh- We obtain 



< \ -= I | 1 ± & | (90) 



where, 



£h — e°h+ir ~ e h + { ( t '>l+7r,/l'+7r,>l+7r,/l'+7r + Vh,h',h,h' 



\ 2 

— Vh+ir,h',h+Tr,h' ~ Vh,h'+ir,h,h'+"!r )V h i 



+ v h+Tr,h',h+TT,h' — Vhh'h.h' 



Aft — —2 (Vh+n,h'+TT,h,h' — Vh+TT,h',h'+TT,h) Uh'Vh' (21) 



With Eq. (120l and Eq. flzT]) we get the self-consistent 'gap' equation 

Aft = —2 (Vh+n,h'-Hr,h,h> ~ Vh+-ir,h',h'+n,h) 7=^= = „ = 

4# ^ Aft/ 

= iV ^ TP 4- A 2 = ' ^ 

Solving the gap equation numerically allows us to calculate the new HF energy as a 
function of g. For iV — > oo, Eq. (1221) gives a nontrivial solution even for g — > 0. 

As a last step we calculate the RPA in the new basis. For this we introduce in analogy 
to Eq. ([9]) an RPA excitation operator of the following form 

Qk = E *>H - Q=\p-h\ (23) 

ph 

Proceeding in a similar way to what we have done in the spherical basis we can find the 
amplitudes X, Y from corresponding RPA equations. We show the spectrum in Fig. El We 
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Figure 5: RPA excitation spectrum in deformed basis in the case of N = 180. 

see that the spectrum now contains a gap corresponding to Eq. ( !2T|) and is qualitatively 
similar to the exact solutions for an anisotropic Heisenberg Hamiltonian (g ^ 1) [3]. 

The fact that the spectrum in Fig. [5] does not resemble much the exact one is not 
completely surprising, since we know that it is built on qualitatively wrong (symmetry 
broken) groundstate. However, the relatively good performance of the HF-RPA scheme 
in the symmetry unbroken phase motivates us to go beyond the standard RPA approach 
in applying symmetry restoration techniques pj] and Self-Consistent RPA (SCRPA) [5]. 
The latter has already produced an interesting result for the 1D-AFH model in Ref. [6]. 

In this work we made a preliminary application of SCRPA to some cases of finite 
number of sites. Without going into details, let us shortly repeat the principles of SCRPA. 
It essentially consists in evaluating Eq. ffTUl) not with the HF ground state as in standard 
RPA but with a ground state containing RPA correlations, as this is the objective proper 
of RPA. In this exploratory application to the AFH model, we found it most convenient 
to solve the vacuum condition Q u \0) = 0, mentioned above, pertubatively. Neglecting 
3p-3h and higher configurations one obtains [IB] : 



Evaluating Eq. (1101) with Eq. (I24p yields an RPA matrix which depends in a non-linear 




(24) 



with 



■p\h\,p2h,2 [ \Plh\,p2h2 




(25) 
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Figure 6: Excitation spectrum u(q) as function of the momentum transfer, for N = 6. 
Exact results are indicated by black squares, RPA results by blue stars and SCRPA results 
by red circles. 



way on the amplitudes X,Y. Solving this non-linear system of equations is numerically 
not completely trivial and consequently in a first attempt we restricted ourselves to quite 
small system sizes. In Fig|H]we compare the results from an exact diagonalization for the 
six sites problem with half filling with the those from SCRPA. We see that there is close 
agreement and a strong improvement of SCRPA over standard RPA. In Fig. [7] we show 
the same but as function of the coupling constant g. Again we see the very important 
improvement of the SCRPA specially around the phase transition point g — 1. In Fig[8]we 
also show the SCRPA results for the case of 20-sites at half filling. We see that the lowest 
branch in SCRPA is strongly lifted up with respect to standard RPA. Though we do not 
have the exact solution at hand for this case, our results may indicate that in SCRPA no 
artificial gap is present any longer in the thermodynamic limit. In order to confirm this 
we have to solve SCRPA in the thermodynamic limit or, at least, approaching it. This 
shall be an investigation for the near future. 

In summary, we evaluated for the first time the full RPA solution for the ID Anti- 
Ferromagnetic Heisenberg model in the Jordan- Winger representation. It is shown that 
one obtains interesting results for spectral functions and sum rules in the symmetry un- 
broken phase in spite of the fact that a low-lying discrete state gets (slightly) detached 
from the continuum what is unphysical. The not translationally invariant Hartree-Fock 
solution still lowers the energy but the corresponding RPA shows, not unexpectedly, a 
strong artificial gap in the spectrum. The encouraging results of standard RPA motivated 
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Figure 7: Excitation spectrum u(q) as function of the anisotropy g, for N = 6 and for 
all momentum transfer. Exact results are indicated by black solid line, RPA results by 
blue dotted line and SCRPA results by red dashed line. For g > 1 the RPA results are 
obtained in the deformed basis. 

us to make a first application of an improved version of RPA, the so-called Self-Consistent 
RPA (SCRPA) [3] for the time being only for small systems. Strong improvement of the 
results of SCRPA over standard RPA can be observed. One may speculate that SCRPA 
cures the artificial gap problem seen in Fig. [3] for the structure function in the thermody- 
namic limit. The latter is, however, a non trivial numerical problem in SCRPA which we 
will try to solve in future work. 
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Figure 8: Excitation spectrum u(q) as function of the momentum transfer for N = 20 
(also for the thermodynamic limit). For the exact result, in the thermodynamic limit, we 
present the lower and upper boundary of the spectrum, they are indicated in black color 
(solid line). For the RPA results, in the thermodynamic limit, we present the collective 
mode and both boundaries of the continuum states; they are indicated in black color 
(dashed line). In the discrete case, RPA results are indicated by blue stars and SCRPA 
results by red circles. 

References 

[1] W. Heisenberg, Z. Phys. 49, 619 (1928). 
[2] H. Bethe, Z. Phys. 71, 205 (1931). 

[3] H.-J. Mikeska and A. K. Kolezhuk, Lecture Notes in Physics vol.645 (Berlin: 
Springer), p. 1-83 (2004). 

[4] A. H. Bougourzi, M. Couture and M. Kacir, Phys. Rev. B 54, R12669 (1996); M. 
Karbach, G. Miiller, A. H. Bougourzi, A. Fledderjohann, and K.-H. Mutter, Phys. 
Rev. B 55, 12 510 (1997). 

[5] J.G. Hirsch, A. Mariano, J. Dukelsky, P. Schuck, Ann. Phys. 296, 187 (2002); A. 
Rabhi, R. Bennaceur, G. Chanfray, and P. Schuck, Phys. Rev. C 66, 064315 (2002); 
A. Storozhenko, P. Schuck, J. Dukelsky, G. Ropke, A. Vdovin, Ann. Phys. 307, 308 
(2003); M. Jemai, P. Schuck, J. Dukelsky and R. Bennaceur, Phys. Rev. B 71, 085115 
(2005). 



12 



[6] P. Kriiger and P. Schuck, EuroPhys. Lett. 72, 395 (1994). 

[7] M. Takahashi, Phys. Rev. B 40, 2494 (1989). 

[8] P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

[9] S. Rodriguez, Phys. Rev. 116, 14741477 (1959); E. Lieb, T. Schultz, and D. Maltis, 
Ann. Phys. (N.Y.) 16, 407 (1961); L. N. Bulaevski, Sov. Phys. JETP 16, 685 (1963). 

[10] D.J. Rowe, Phys. Rev. 175, 1283 (1968); Rev. Mod. Phys. 40, 153 (1968). 

[11] P. Ring and P. Schuck, The Nuclear Many-Body Problem, Springer, Berlin (1980) 

[12] D. J. Thouless, Nucl. Phys. 22, 78 (1961). 

[13] T. Todani and K. Kawasaki, Prog. Theor. Phys. 50, 1216 (1973). 

[14] G. Muller, H. Thomas, H. Beck, and J. C. Bonner, Phys. Rev. B 24, 1429 (1981). 

[15] T. Schneider, E. Stoll, and U. Glaus, Phys. Rev. B 26, 1321 (1982). 

[16] H. Beck and G. Muller, Solid State Commun. 43, 399 (1982). 

[17] A. Auerbach, Interacting Electrons and Quantum Magnetism, Springer- Verlag, New 
York (1994). 

[18] D. S. Delion, P. Schuck and J. Dukelsky, Phys.Rev. C 72, 064305 (2005); D. Gam- 
bacurta, M. Grasso, F. Catara, and M. Sambataro, Phys. Rev. C 73, 024319 (2006). 



13 



